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The paper contains a rigorous proof of the absence of quasi-long-range order in the random-field 
O(N) model for strong disorder in the space of an arbitrary dimensionality. This result implies 
that quasi-long-range order inherent to the Bragg glass phase of the vortex system in disordered 
superconductors is absent as the disorder or external magnetic field is strong. 
The nature of the vortex phases of disordered superconductors is a subject of active current investigations. A 
plausible picture includes three phases Q : vortex liquid (VL) , vortex glass (VG) and Bragg glass (BG) . In all those 
phases an Abrikosov lattice is absent Q and only short-range order (SRO) is expected in VG and VL. A higher degree 
P-h ' of ordering is predicted in BG. It is argued that in this state the vortex array is quasi-long-range-ordered ||||. Thus, 
qq , Bragg peaks can be observed in BG (||| as if the system had an Abrikosov lattice. In the other phases Bragg peaks 
are not found ||. The phase transitions from BG to VG and VL are presumably associated with topological defects 

The above picture is supported by variational and renormalization group calculations for the random-field XY 
I | model dJlfiH wn ' cn * s ^ ne simplest model of the vortex array in disordered superconductors. Besides, this model is 
useful for our understanding of many other disordered systems Hence, the detailed knowledge of its properties 
is important. Unfortunately, the only rigorous result about the random-field XY model is the absence of long-range 
! order (LRO) The present paper contains a new rigorous result: QLRO is absent as the disorder is sufficiently 
strong. It is interesting that our proof is quite simple in contrast to the very nontrivial demonstration of the absence 
C °f LRO for arbitrarily weak disorder j| . 

It is expected that the vortex system of the disordered superconductor has no ordering for strong disorder due to 
appearance of the dislocations. The XY model is a convenient polygon for understanding of the role of the topological 
defects. Since the vortex glass phase of disordered superconductors corresponds to strong disorder, our rigorous result 
q ' supports the conjecture that the vortex glass phase has no QLRO. 

The simplest model of the vortex array in disordered superconductors has the following Hamiltonian 

: , 

>■ H = / d 3 r[K(Vu(r)) 2 + hcos(2%u(r)/a- 6(r))], (1) 

0\ ; J 

I ' where u is the vortex displacement, a the constant of the Abrikosov lattice in the absence of disorder, the random 
phase. The one-component displacement field u(r) describes anisotropic superconductors. The generalization for the 
' isotropic case is straightforward. The ordering can be characterized in terms of the form-factor 

o : 

G(r) = (cos27t(m(0) - u(r))/a), (2) 

where the angular brackets denote the thermal and disorder average. This correlation function contains information 
£^ about neutron scattering. LRO corresponds to a finite large-distance asymptotics G(r — > oo) — ► constant, QLRO is 
described by the power law G(r) ~ r~ v and SRO corresponds to the exponential decay of the correlation function 
G(r) at large r. 

We demonstrate that as the random-field amplitude h Eq. (Q) is large the system possesses only SRO. Since h 
depends on the strength of the disorder in the sample and the external magnetic field B we see that SRO corresponds 
to the situation at which either the disorder or magnetic field is strong. We consider not only the XY model (|l|) but 
• j—{ , also the more general random-field O(N) model. Its Hamiltonian has the following structure 

H = -jVJSiSj - H^SiHi, (3) 

m 

where Si are the p-component unit spin vectors on the simple cubic lattice in the D-dimensional space, rii is the 
random unit vector describing orientation of the random field at site i, the angular brackets denote the summation 
over the nearest neighbors on the lattice. The XY Hamiltonian corresponds to p = 2. In this case the relation 
between (@) and © is given by the formulae S x = cos27rit/a, S y = sin27ru/a, where S x ,y are the spin components. 

The idea of our proof is based on the fact that the orientation of any spin Si depends mostly on the random fields 
at the nearest sites as the random-field amplitude H is sufficiently strong. We shall see that the knowledge of the 
random fields in the region of size Nb with the center in site i, where b is the lattice constant, allows us to determine 
the orientation of the spin Si with the accuracy exp(— constantV). Thus, the orientations of any two distant spins 
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depend on the realizations of the random field in two non-intersecting regions up to exponentially small corrections. 
The values of the random field in these regions are uncorrelated. Hence, the correlations of the distant spins are 
exponentially small. 

Below we consider the case of the zero temperature. Hence, the system is in the ground state. We assume that the 
amplitude of the random field 

H = 2DJ{1 + V2 + 5),S> 0. (4) 

Let us estimate the angle between an arbitrary spin Si and the local random field hi = ifiij. The Weiss field 
Hw = hj + Hj acting on the spin includes the random field hi and the exchange contribution Hj = JY]j Sj, where 
Sj are the nearest neighbors of the spin Sj. The latter Hj < 2D J since the number of the nearest neighbors is 2D, 
where D is the spatial dimension. Hence, the minimal possible modulus of the Weiss field is (H — Hj) > (H — 2D J). 
Any spin is oriented along the local Weiss field %. Let us consider the triangle two sides of which are hi and H / , 
and the third side is parallel to Si . The laws of sinuses allow us to show that the maximal possible angle between h; 
and Si is 4>\ — &rcsm(2DJ/H). 

We shall now determine the orientation of the spin Si iteratively. Let the zero approximation S? be oriented along 
the random field hi. Let the first approximation Sj be oriented along the Weiss field H^,-, calculated with the zero 
approximation for the neighboring spins: H-%r = hi + J^2j n?, where • denotes the summation over the nearest 
neighbors. The second approximation Sf is determined with the Weiss field in the first approximation, etc. Any 
approximation Sf depends on the random fields only at a finite set of the lattice sites. The distance between any 
such site and site i is no more than kb, where b is the lattice constant. 

In any approximation the Weiss field > H — 2D J. Let be the difference between the fcth and (k — l)th 

approximations for the spin Si. Then |H{^ — Hj^ < 2DJm k , where m k denotes the maximal value of sf. Hence, 
we find with the laws of sinuses that the angle between S k and s\ k ^ is less than (f>k = arcsin(2Z? Jm k /{H — 2D J)). 
Since sf +1) < 2sin(0 fe /2), one obtains the following estimation: 

m (fc+1) < ^/ 2{1 - [2DJm k /(H -2DJ)} 2 } = ^2[2DJm k /(H-2DJ)] 2 /[l + yd" [2DJm k /(H - 2D J)) 2 } < 

2V2DJm k /{H - 2D J) = m k /[l + S/V2] < m 1 /[l + S/V2] k , (5) 

where Eq. (||) is used. 

Now we are in the position to estimate the correlation function (j^). In terms of the 0(N) model it is given by the 
expression 

G(r) = (S(O)S(r)). (6) 

Let N — [r/2b] — I, where the square brackets denote the integer part. We decompose the values of the spins in the 
following way S(x) = S jv (x) + J2k>N sk { x )- The iVth approximations S Ar (x) and 8^(0) depend on the orientations 
of the random fields in different regions which have no intersection. Hence, the correlation function (S Af (0)S JV (r)) is 
the product of the averages of the two multipliers S N and thus equals to zero due to the isotropy of the distribution 
of the random field. All other contributions to Eq. (^|) can be estimated with Eq. (||) and are exponentially small as 
the functions of N. This proves that the correlation function G(r) has an exponentially small asymptotics at large r. 
Thus, for strong disorder (||) both LRO and QLRO are absent. 

A challenging question concerns the presence of QLRO in the weakly disordered random-field systems. Another 
interesting related problem is the question about QLRO in the random-anisotropy O(N) model [[10|. Unfortunately, 
our approach cannot be directly generalized for this problem since the zeroth-order approximation S? is not unique 
in the random-anisotropy model: any spin has two preferable orientations. 

In conclusion, we have proved that for strong disorder the random-field O(N) model has no QLRO. The renormal- 
ization group calculations jl0| suggest that at N > 2 QLRO is absent for arbitrarily weak disorder. On the other hand, 
in the random- field 0(2) model without dislocations the renormalization group |I| predicts QLRO. Our result shows 
that in the system with the topological defects QLRO is absent at least as the disorder is strong. This prediction is 
relevant for the vortex glass state of the disordered superconductors. 
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